The structure of the central disk of NGC 1068: a clumpy disk model 
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ON . NGC 1068 is one of the best studied Seyfert II galaxies, for which the blackhole 



mass has been determined from the Doppler velocities of water maser. We show that the 



standard a-disk model of NGC 1068 gives disk mass between the radii of 0.65 pc and 

1.1 pc (the region from which water maser emission is detected) to be about 7xl0 7 Mq 

(for a = 0.1), more than four times the blackhole mass, and a Toomre Q-parameter for 
d ■ 

the disk is ~0.001. This disk is therefore highly self-gravitating and is subject to large- 



amplitude density fluctuations. We conclude that the standard a-viscosity description for 



H ! the structure of the accretion disk is invalid for NGC 1068. 

cz i 

In this paper we develop a new model for the accretion disk. The disk is considered 
to be composed of gravitationally bound clumps; accretion in this clumped disk model 
arises because of gravitational interaction of clumps with each other and the dynamical 
frictional drag exerted on clumps from the stars in the central region of the galaxy. The 
clumped disk model provides a self-consistent description of the observations of NGC 1068. 
The computed temperature and density are within the allowed parameter range for water 
maser emission, and the rotational velocity in the disk falls off as r -035 . 
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1. Introduction 

The discovery of emission from water masers in the central region of Seyfert II galaxy 
NGC 1068 provides an opportunity to study the properties of the associated accretion disk. 
The mass of the central black hole is estimated to be 1.5xl0 7 M© (Greenhill & Gwinn, 
1997) from the measurement of the Doppler velocities of the masing spots. 

The nucleus of NGC 1068 is heavily obscured, and the luminosity of the central source 
is determined from the observed flux by modeling the dust obscuration and the scattering 
of photons by ionized gas in the nuclear region. According to a careful analysis carried 
out by Pier et al. (1994) the bolometric luminosity of NGC 1068 is estimated to be about 
8xl0 44 erg s _1 ; the luminosity is perhaps uncertain by a factor of a few. 

We use the blackhole mass and the bolometric luminosity to construct the standard 
viscous disk model for NGC 1068 (§2) and find that the disk is highly self-gravitating 
thereby rendering the a-disk model inapplicable. The effect of the irradiation of the disk 
from the central source does not modify this conclusion (§2b). 

In section 3 we present a new model for the disk in NGC 1068 composed of gas clumps. 
The accretion in this case arises as a result of gravitational interaction amongst clumps 
(§3.1); the structure of the clumped-disk model (CD model) of NGC 1068 is described in 
§3.2. 

The velocity of the masing spots appears to be falling off with distance from the center 
as r -0 ' 35 (Greenhill & Gwinn, 1997), which is less rapid than the Keplerian power law of 
—0.5. One possible reason for this could be that the disk of NGC 1068 is sufficiently 
massive so as to modify the rotation curve. However, we show in §3 that this is not so. 
Another possibility is that the flattening of the rotation curve is due to a central star 
cluster. The influence of the star cluster on the accretion rate and the disk structure is 
discussed in section 3.2. 

2. Standard thin disk model for NGC 1068 

The standard viscous accretion disk model for NGC 1068, including the irradiation of 
the disk from the central source, is presented below. Throughout this paper we take the 
mass of the blackhole at the center of NGC 1068 to be 1.5xl0 7 M© (Greenhill & Gwinn, 
1997), and the bolometric luminosity to be 8x10 erg s _1 (Pier et al., 1994). 

The theory of thin accretion disk is well developed and is described in a number 
of review articles and monographs e.g. Frank, King & Raine (1992). When the flux 
intercepted by disk is not small compared to the local energy generation rate, such as 



that expected for the masing disk of NGC 1068, the incident flux must be included in 
determining the thermal structure of the disk. 

The fractional luminosity intercepted by the disk depends on the disk geometry, the 
scattering of radiation by the coronal gas etc. For instance if the dominant source of 
radiation intercepted by the disk were the scattered radiation from an extended corona 
then we might expect the incident flux at the disk to be roughly uniform. On the other 
hand for a flaring or a warped disk the radiation intercepted from the source directly 
might dominate, and the incident radiation in this case depends on the inclination angle 
of the normal to the disk. We adopt the second model to analyze the disk structure; it is 
straightforward to consider other possibilities, but we do not pursue these since the main 
result of this section viz., the standard a-disk model for NGC 1068 is highly unstable to 
gravitational perturbations turns out to be independent of whether we include irradiation 
and radiative forces in the calculation of disk structure. 

Let the luminosity of the central source be L c . The flux incident at the disk, i^ n , a 
distance r from the central source is taken to be 

p. ( ) - ^ c ^H _ /3HL C 
4 7rr 2 ^ r 47rr 3 

where H is the vertical scale height, the constant factor (3 is defined by dH/dr = (3H/r, 
and determines the fraction of the flux that is intercepted by the disk. The upward energy 
flux due to mass accretion rate M is 

F up (r) = ^ 2 M, (2) 

and the ratio of these fluxes Fi n /F up rs e[3H/R sc , where e ~ 0.1 is the efficiency of the 
conversion of the rest mass to energy by the central source, and R sc is the Schwarzschild 
radius of the blackhole, O = (GMt/r 3 ) 1 ' 2 is the angular rotation speed, and M t is the 
mass contained inside the radius r. The effective temperature of the disk at r is given by 

aT^ ff (r) = F up + F in = F t . (3) 

The temperature at the mid-plane of the disk (T c ) can be calculated by considering 
the first moment of the radiative transfer equation, and making use of the Eddington 
approximation, and is given by 
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where To = kE is the optical depth of the disk. The opacity k, for masing disks, with 
Hi density and temperature in the range 10 8 -10 10 cm" 3 and 100-1000 K respectively is 
dominated by metal grains and is given by (Bell & Lin 1994) 



K 



0.1T 1/2 cm 2 g _1 . (5) 



The temperature in the disk mid-plane is not sensitive to the height where the energy 
is deposited so long as the optical depth of the disk to the incident radiation is much 
greater than one. 

The solution to the hydrostatic equilibrium equation in the vertical direction, which 
includes the radiation pressure at the disk surface ~ (2Fi n + F up )/c, yields 



p c ci « VZH + ^f[l~ ^f ) ■ (6) 

where E is the surface mass density, c s & p c are the sound speed and gas density in the mid 
plane of the disk. Taking the effective viscosity in the disk to be ac s H the mass accretion 
rate M is given by 

27raEc 2 . ._. 

M = -jj-l. (7) 

Equations (4), (6) and (7) are three equations in three unknowns viz. E, H and T c , 
which are solved numerically and the results are shown in fig. 1. Also shown in fig. 1, for 
comparison, are the solutions for (3 = i.e. when the irradience of the disk is neglected. 
The solution in the latter case can be obtained analytically and is given below 

T c = 200 a-^m^M^r- 1 K, (8) 

n = l^a-^m^M^r-^ 2 cm" 3 , (9) 

and the Toomre-Q parameter for the stability of the disk is 

Q = -% = 2.4x lO-V^m-^M^V 3 / 2 , (10) 

where m is the accretion rate in terms of the Eddington rate, Mi = M/(10 7 Mq), and 
r is the radial distance from the center in parsecs. For the NGC 1068 system Mf rs 1.5 
and m « 0.4. The inner and the outer radii of the masing disk are at 0.65pc and l.lpc 
respectively. 



For the numerical results shown in fig. 1 we took the value of (3 such that the disk 
intercepts about 50% of the flux from the central source. In this case the disk temperature 
is dominated by the incident flux, and the temperature in the disk midplane is close to 
T e ff pa 511K at r = lpc. The structure of the disk (E, Q, H, T c ) in this case is almost 
independent of the opacity of the gas and consequently it is unaffected by any uncertainty 
in k. The molecular mass of the masing disk, when irradiation is included, is ~ 7.OxlO 7 M0 
(for a = 0.1), and the Q is about 10~ 3 (see fig. 1). For the irradiation dominated disk, 
the disk mass decreases with a as a -09 , and the Q increases as a - 9 . Since Q should be 
greater than one for stability, we see that the masing disk of NGC 1068 is highly unstable 
to gravitational perturbation. A decrease in the irradiation flux makes the disk more 
unstable. 

The main conclusion of this section is that the standard a-disk model for the NGC 
1068 masing disk is inconsistent. The disk according to this model is highly unstable, and 
should fragment into clumps. A consistent analysis in this case should include the effect 
of self-gravity, and clump interaction, which is described below. 

3. A model for clumpy self gravitating disks 

A disk with small value for (1 — Q) is likely to develop spiral structure which can 
transport angular momentum outward. However, when Q <C 1, as in the case of NGC 1068 
(see §2), the disk is likely to breakup into clumps, and the accretion rate is determined 
by the gravitational interaction among these clumps. Accretion in self gravitating disks 
was considered by Paczynski (1978), Lin & Pringle (1987), Shlosman & Begelman (1987), 
Shlosman et al. (1990), and has been investigated more recently by Kumar (1998) in some 
detail for clumpy disks. In §3.1 we derive the results that we need to construct model for 
clumpy disks (CDs), and its application to NGC 1068 is discussed in §3.2. 

§3.1 Velocity dispersion and accretion rate in a clumpy disk 

Consider a disk consisting of clumps of size / c , and take the mean separation between 
clumps to be d c . The tidal radius of a cloud, the distance to which the gravity of the cloud 
dominates over the gravity of the central mass, is dt ~ r(m c /M t ) 1 ' 3 ; where m c ~ a c l 2 c is the 
cloud mass, a c is the surface mass density of the cloud and Mt is the total axisymmetrically 
distributed mass contained inside the radius r. 

The change to the displacement amplitude of the epicyclic oscillation of a clump as 
a result of gravitational interaction with another clump, with impact parameter d > dt, 
can be shown to be ~ m c r 3 /(Mtd 2 ). Gravitational encounters with d < dt are almost 
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adiabatic and these leave the epicyclic energy of clumps unchanged (two particles moving 
on circular orbits merely interchange their trajectories in such an encounter). Thus the 
dominant gravitational interactions for exciting epicyclic oscillation are those with impact 
parameter of ~ dt, and the change to the epicyclic amplitude in such an encounter is 

5r ~ m c r 3 / (M t df) ~ dt, or 

^ /v/3 (^) 1/3 - <"> 

where 

M r = ^^l=7rr 2 a(r), (12) 

and a(r) is the mean surface mass density of the disk at radius r. 

The mean time interval for a clump to undergoes strong gravitational interaction with 
another clump i.e. with impact parameter dt is ti nt ~ Vt~ 1 {d c /dt) 2 , and so long as the cloud 
size (l c ) is not much smaller than d c the timescale for physical collision between clumps 
is also of the same order i.e. clouds undergo physical collision after having undergone 
one strong gravitational encounter on the average, and so the cloud velocity dispersion is 
~ QSr or 



(rfi) 



M r 



1 1/3 r 4 12/3 



(13) 



nM t _ 

The expression for v r is similar to that given by Gammie et al. (1991) for the velocity 
dispersion of molecular clouds in the Galaxy. 

We assume that the kinetic energy of the relative motion of colliding clouds is dis- 
sipated and their orbit is circularized. Colliding clumps might coalesce so long as the 
size of the resulting object does not exceed the maximum length for gravitational insta- 
bility, Imax ~ (rM r )/M t ; clouds of larger size are susceptible to fragmentation due to the 
rotational shear of the disk which limits their size. 

The characteristic time for a clump to fall to the center can be estimated using equation 
(11) and is given by 

■—•o'er 

and the associated average mass accretion rate is 



MrV' 3 Ur^' 3 



M,n(r)M r ^j ^j . (15) 

This result is applicable so long as the density contrast between the clump and the inter- 
clump medium is about a factor of two or larger, and l c is of order d c . 

For l c <C d c clouds undergo several gravitational encounters with other clouds before 
undergoing a physical collision, and during this time the random velocity of clouds con- 
tinues to increase. The rate of increase of velocity dispersion in two body gravitational 
interaction is given by 

dv 2 _ G 2 m c M r n 

dt j-2^2 

We assume that the scattering gives rise to isotropic velocity dispersion, and so the scale 
height for the vertical distribution of clumps is ~ v/Q. The time for a clump to undergo 
physical collision, assumed to be completely inelastic, is t co \ ~ fi _1 (<i c /Z c ) 2 , and therefore 
the velocity dispersion is 

v 2 ~ ^1. (17) 

This corresponds to the Safronov number being equal to one. Since the cloud collision time 
is much greater than O -1 the effective viscosity is suppressed by a factor (Qt co i) 2 compared 
to the case where collision frequency is greater than O, (cf. Goldreich and Tremaine, 1978), 
i.e. v e ~ v 2 (t~ ol /VL 2 ), and thus the mass accretion rate is given by 

M A fk\ ~,, /M r \ 4/3 fd c \ 2/3 (l c 



This accretion rate is smaller than given by equation (15) by a factor of (dt/l c ). 

We can parameterize the effect of unknown size distribution and separation between 
clumps on the accretion rate and the velocity dispersion by a dimensionless parameter r\ 1 
and rewrite equations for M and v r in the following form for future use: 

M = T]n(r)M r (r) ( -^- j , (19) 



K nM t 
and 



<Y = nrn(r) ( -^ ) . (20) 



For l c ~ d c ~ Imaxi we see from equations (16) and (18) that rj ~ 1. The effective 
kinematic viscosity in this case is approximately /^ ax fi ~ Q~ 2 H 2 Vt (where Q ~ l max /H is 
the Toomre Q-parameter, i7 ~ v r /fl), same as in the ansatz suggested by Lin & Pringle 
(1987) for self-gravitating-disk. We note that the relative velocity of collision between 
clouds is smaller than their orbital speed by a factor of the ratio of the total mass to the 
molecular mass, and as long as the cooling time for post-collision gas is less than the time 
between collisions the clouds are not smeared away due to heating and the disk remains 
clumpy. 

§3.2 Application to NGC 1068 

The equation (19) can be recast in the following form 

M 2 

'""'■ -4>, (21) 



M r + M c 

where 

7T 2 rM 2 /3 , v 

The solution to this equation in the case where the central mass dominates is 

M r * V 1/2 (M C V 2 + tl^ , (23) 

and from equation (20) we find the velocity dispersion of clouds to be 

v r ^r] 2/3 [GM] 1/s . (24) 

For NGC 1068, the luminosity ~ 8xl0 44 erg s _1 , M c ~ 1.5xl0 7 M , and M ~ 8xl0 24 
g s _1 . Substituting these numbers in the above equations we find 

M r w 1.3 x 10 6 M Q r 1/2 ?7" 1/3 , (25) 

and 

t) r «8kms-S 2/3 , (26) 

where r is measured in parsecs. Note that the disk mass in this model is more than an 
order of magnitude smaller than in the standard ct-disk model discussed in §2. The velocity 
dispersion of clumps is larger than the sound speed, and is of order the observed scatter 
in the velocity of masing spots in NGC 1068. 
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The number density of molecules and the gas temperature can be obtained by solving 
the hydrostatic and the thermal equilibrium equations for blobs in the vertical direction 
in the presence of incident flux (see §2). We find the thermal temperature of blobs to be 
about 510 K at r = lpc (see fig. 1), and the density scale height H ~ 1.5 x 10 16 cm. The 
mean number density of H2 molecules n ~ M r /(47Tr 2 H[i) = 5 x 10 8 cm -3 . Thus both 
the number density and the temperature are within the allowed range for water maser 
emission. 1 The size of clumps is < r(M r /M t ) ~ 0.1 pc and their mass is ~ 10 3 M . Since 
the Jeans mass is ~ 10 M the clumps could have some star formation activity. However, 
the efficiency of star formation is usually quite low, of order a few percent even for the 
giant molecular clouds, so we don't expect the gas in the clumps to be converted into stars 
in the short lifetime of order 10 6 years or less for the clumps in the masing disk of NGC 
1068. 

Note that the velocity dispersion of blobs is independent of r (eq. [26]), so the scale 
height for the vertical distribution of blobs increases as r 3 / 2 provided that the velocity 
distribution of clouds is nearly isotropic; this increase of disk thickness with r is more 
rapid than in the standard ct-disk model considered in the last section. At a distance of 2 
pc from the center of NGC 1068 the scale height is about 0.3 pc (for rj = 5), so the flaring 
disk blocks a significant fraction of the radiation from the central source. 

Clouds colliding at a relative speed of 8 km s _1 raise the gas temperature to approxi- 
mately 10 4 K which is too hot for maser emission. However, the cooling time of gas at the 
density of ~ 10 9 cm -3 , is of order 100 years, which is short compared to both the orbital 
time and the collision time of order 10 4 years 2 ensuring that the disk remains cold and 
that a steady state solution exists. 

The slope of the rotation curve in the disk follows from the use of equation (25) 

^ = -0. 8 + I^ = -0. 8 + ^l K -0.5+^. (27) 

d\nr 2 dlnr 4M C 1/2 4M C v ' 

The disk mass for the NGC 1068 system, in our model, is small (M r /M t ss 0.1), and so 
the slope of the rotation curve is very close to —0.5. 

A slower fall off of the rotation curve requires a more rapid increase of the total mass 
(Mt) with r than in the model discussed above. This could arise for instance if there is a 

1 The allowed range for number density for water maser emission is 10 8 — 10 11 cm -3 

and the temperature range is 200 — 1000 K. 

2 The density enhancement of the shocked gas depends on the strength of magnetic field 

in the clumps and is order unity for equipartition magnetic field. 
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star cluster at the center. The rotation curve falls off as r -035 when the mass of the star 
cluster within a parsec of the center of NGC 1068 is about 8xl0 6 M Q . Thatte et al. (1997) 
have carried out near infrared speckle imaging of the central 1" of NGC 1068, and conclude 
that about 6% of the near-infrared light from this region is contributed by a star cluster. 
They estimate that the stellar mass within 1" (~ 50pc) of the nucleus of NGC 1068 is 
about 6xl0 8 Mq; if the density in this cluster falls off as r -2 , as in a singular isothermal 
sphere, then the expected stellar mass inside 1 pc is 1.2xl0 7 M . We note that the stellar 
mass within 100 pc of our own galactic center is estimated to be in excess of 5xl0 8 M Q 
(Kormendy & Richstone 1995, Genzel et al. 1994), and the mass enclosed within radius 
r is seen to increase almost linearly with r for r > 2pc. Thus the possibility of a similar 
stellar mass cluster at the center of NGC 1068 is not surprising. 

We discuss below the effect a star cluster has on the structure of the accretion disk. 

§3.2.1 Effect of a star cluster on accretion disk 

Let us consider that the stellar mass within a radius r of the center is M s (r). The disk 
mass, as before, is taken to be M r (r), and the central mass is M c . The use of equation 
(23), which still applies with M c replaced by (M c + M s ), implies that the disk structure 
is not much affected at small radii where M s (r) <C M c . At larger radii, where the stellar 
mass becomes comparable to or exceeds M c , the mass of the gas disk (M r ) increases with 
r as r x ' 2 [M c + Ms) 1 ' 2 , which is somewhat more rapid than the case of M s = considered 
in §3.1. However, for the masing disk of NGC 1068 M s < M c for r < 1.6 pc, and the effect 
of the stellar cluster on the disk structure, in particular the number density of molecules 
and the gas temperature, is small. 

The velocity dispersion of clouds due to gravitational encounters is proportional to 
the local surface mass density a of gas and is unaffected by the stellar cluster (see eq. [24]). 
Therefore, the disk thickness (H) at first increases with radius as r 3 ' 2 and then flattens 
out when M s starts to dominate the total mass: 

H ^J^_ „ Vr 3 a(r) 

H ~n(r)~ M t ■ {2H) 

For a oc r -3 ' 2 , expected of the CD model with constant M, H/r ex r 1 ' 2 , however a more 
rapid decrease of a with r leads to a corresponding decrease of H/r. 3 



3 A decrease of a with r that is more rapid than r 3 / 2 leads to a drop in M. How- 
ever, the assumption of steady state accretion breaks down beyond some radius where the 
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The accretion rate calculated in §3.1 is modified due to the dynamical friction suffered 
by clumps in the disk by the star cluster; the orbits of stars in the cluster are gravitationally 
perturbed by the clumps so that on the average the density of stars behind a clump is 
greater than the density in the front. We estimate below the accretion rate that arises as 
a result of the frictional drag exerted by stars. We assume that clumps are not stretched 
out by the tidal force of stars, which is valid so long as the cloud mass density is greater 
than the mean mass density associated with the star cluster i.e. n# 2 > M s /(mH 2 r3 ) or 
n H2 >2x 10 8 cm" 3 (M s /10 7 M Q )(r/lpc)- 3 ; this condition is satisfied for NGC 1068. 

The dynamical friction timescale for a clump to fall to the center is (cf. Binney and 
Tremaine, 1987) 

r 2 M t v rot 
df Gm c M s In A' [U) 

where m c is clump mass, A ps rv 2 ot /(Gm c ), and v 2 ot = GM t /r. Taking the clump size to 
be the largest length scale for gravitational instability in a shearing disk i.e. l c ~ n 2 Ga/Q 2 , 
we find m c « 7rM 3 /M t 2 . Thus, In A w 31n(M t /M r ) is about 7 for the NGC 1068 system. 
Substituting these in the above equation we find 

'* * To" \j^mI) ■ (30) 

For M s « M t /S and M r /M t « 0.1, values applicable to the NGC 1068 system, the dy- 
namical friction time is of the same order as the timescale for clumps to fall to the center 
due to gravitational encounter with other clumps (see eq. [22]). Thus we see that a star 
cluster of modest mass at the center of NGC 1068 can both give rise to a sub-Keplerian 
rotation curve, as perhaps observed by the water maser, and also because of its dynamical 
friction on clumps remove angular momentum of the gaseous disk resulting in accretion 
rate that is of the same order as needed to account for the observed luminosity (the total 
M is the sum of the accretion rate due to gravitational interaction between clumps, and 
the dynamical friction drag exerted by stars). 

It was pointed out by Ostriker (1983) that even a smooth disk is subject to frictional 
drag from a star cluster. He showed that the characteristic drag time for removing angular 
momentum of a gaseous disk is of the same order as the relaxation time for the star cluster. 
The accretion rate due to this process, in a form applicable to the NGC 1068 system, is 
given by 



accretion time, O 1 (M C + M s ) 2 /M 2 , is longer compared with the evolution time scale of 
the disk. 
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M~M r Q/ (y) 2 (— ), (31) 

where r* is the radius and m* is the mass of stars in the cluster, and Jo ~ 10 is a 
dimensionless quantity that depends on the ratio of the escape velocity at the surface of a 
star to the stellar velocity dispersion. We see that the accretion rate resulting from friction 
drag from stars on a smooth disk is much smaller than the rate resulting from the frictional 
drag on clumps calculated earlier. 

4. Conclusion 

We find that the structure of masing disk of NGC 1068, as determined using the 
standard a-viscosity prescription, is both inconsistent with the observations and is self 
contradictory. In particular, the disk mass according to this model is 7xl0 7 M , which 
is much greater than the blackhole mass, and the Toomre Q parameter has an extremely 
low value of ~ 10 -3 (for a = 1). This makes the disk highly unstable to gravitational 
perturbations, and suggests that it breaks up into clumps, a conclusion that is inconsistent 
with the smooth disk assumption of standard a-disks. 

We have described in this paper a different model for the disk. We consider the disk 
to consist of gas clumps which undergo gravitational interactions with one another leading 
to an inward accretion of gas. The mass of the masing disk of NGC 1068 according to 
this model is about 10 6 M , and the velocity dispersion of clumps is about 10 km s _1 
which is in rough agreement with the observed velocity dispersion of masing spots. The 
temperature and the density of clumps in this model are approximately 510 K and 5xl0 8 
cm -3 respectively, which are within the allowed parameter range for maser emission. 

However, the rotational velocity of clumps, in this model, are close to the Keplerian 
value, whereas the observations of masing spots indicate a slower fall off of their rotational 
velocity (Greenhill & Gwinn, 1997). One obvious way these results can be reconciled is if 
the mass in stars within 1 pc of the center of NGC 1068 is of order the blackhole mass. 
Speckle observations of the central 1" region of NGC 1068 in near infrared in fact suggests 
that the stellar mass contained within 50 pc of the center is about 5xl0 8 M (Thatte et 
al. 1997), which is sufficient to explain the flattening of the observed rotation curve in the 
masing disk. The dynamical friction exerted by this star cluster on the clumps in the disk 
removes angular momentum at a rate that is of the same order as needed for the nearly 
Eddington accretion rate for the system. Thus the clump-clump gravitational interaction 
and the dynamical friction drag force on clumps together determine the disk structure, 
which we find to be consistent with observations. 

13 



The clumpy-disk (CD) model considered in this paper for the disk of NGC 1068 should 
also apply to any AGN at a distance of about a pc or more from the center, depending 
on the accretion rate and the luminosity, where the disk becomes self-gravitating and the 
standard a model is too inefficient to account for the accretion rate (Shlosman & Begelman, 
1989). 
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Figure Caption 

FIG. 1. — The panel on the left shows the gas temperature (T c ) and the number density of 
B.2 molecules (n c ) in the mid plane of the disk for the standard a-disk model without any 
irradiation from the central source (thin continuous line, and thin dashed line respectively). 
The thick continuous and dashed lines are T c and n c respectively for an a-disk model that 
includes irradiation from the central source. The irradiation is taken to be proportional to 
H/r, where H is the scale height of the disk, and the proportionality constant (/?) is chosen 
so that the flux intercepted by the disk at distance r is about 50% of the flux from the 
central source at r (such a large value for the flux intercepted by the disk corresponds to 
(3 = 50 which might arise if the disk were extremely warped). The panel on the right shows 
the Toomre Q-parameter for the two disk models; thin continuous curve is for the standard 
a-disk with no irradiation, and the thick curve is for a disk that intercepts about 50% of 
the flux. Note that Q is much less than one in both of these models which corresponds to 
the disk being highly unstable to gravitational perturbations. In both of these models we 
chose a = 0.1 (Q scales as roughly a 09 when irradiation dominates), the blackhole mass 
M — 1.5 x 1O 7 M , and the mass accretion rate (m) in units of the Eddington rate was 
taken to be 0.4. 
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